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ABSTRACT 


Gas  hydrate  equilibrium  curves  for  the  ternary  system 
N2-CH4-H2O  were  obtained.  The  hydrate-forming  conditions  of 
pure  nitrogen,  pure  methane  and  of  six  mixtures  of  nitrogen 
and  methane  have  been  investigated  in  the  temperature  interval 
32°F  to  72°F  and  between  pressures  of  300  and  5100  psia. 

The  pressure-composition  equilibrium  diagrams  at  32°F, 
39.7°F  and  44°f  for  the  ternary  systems  were  also  obtained. 

The  information  thus  obtained  was  used  to  evaluate  vapour-solid 
equilibrium  constants  for  methane  and  nitrogen.  With  the  set 
of  equilibrium  constants  so  obtained  curves  representing  the 
equilibrium  constants  for  methane  were  extended  in  the  region 
where  the  values  of  equilibrium  constants  for  methane  are  less 
than  unity. 

The  solid-solution  theory  of  Platteeuw  and  Van  der 

Waals  was  applied  to  predict  the  pressure-composition  diagrams 

o 

for  the  system  at  32  F.  The  theory  was  also  applied  to  predict 
pressure-composition  diagrams  for  the  ternary  systems  CH4-C2H4-H2O 
and  CH^-C2Hg-H20 .  The  phase-diagrams  thus  obtained  also  gave 
the  solid  vapour  equilibrium  constants  for  CH^.  The  experimental 


values  of  the  equilibrium  constant  fell  below  the  theoretical 


values  at  pressures  of  about  100  psia,  coincided  with  the 
theoretical  values  at  about  400  psia,  and  fell  about  25  per 
cent  above  the  theoretical  values  at  pressures  in  the  region 
of  1500  psia. 
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As  a  result  of  this  work  it  is  established  that  gas- 

hydrates  grow  in  two  structures,  referred  to  as  Structure  I  and 

Structure  II.  In  these  structures,  gas  molecules  are  entrapped 

in  cavities  formed  by  water  molecules  linked  together  by  hydrogen 

bonding.  There  are  two  types  of  cavities  in  each  structure  and 

one  of  them  is  considerably  smaller  than  the  other.  Depending 

on  the  size  of  the  molecule  (maximum  molecular  dimension)  the 

gas  will  occupy  the  smaller  or  larger  cavities.  Details  of 

(28) 

the  cavity  geometry  are  given  below  .  Figure  1  and  2  describe 
the  shape  of  the  cavities. 


Smaller  Cavities 

Structure  I 
Structure  II 


( Pentagonal 


dodecahedra) 

Average  diameter 
5.10  A° 

5.00  A° 


Large  Cavities 

Structure  I 
Structure  II 


(Tetrakai  decahedra) 

5.80  A° 
6.70  A° 


Dimensions  of  the  cavities  immediately  restrict 
the  number  of  gases  which  can  form  hydrates.  Any  gas  having 
a  molecular  dimension  larger  than  6.7  A°  cannot  form  hydrates. 
Theoretically  there  cannot  be  any  restriction  on  the  smaller 
molecules.  However,  it  is  found  that  smaller  molecules  like 
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hydrogen  and  helium  are  not  retained  within  the  frame  very 
strongly.  Gases  with  a  hydrogen  atom  in  a  polar -bond  do  not 
form  hydrates.  Halogen  acids,  alcohols,  amines,  organic  acids, 
amides  and  ammonia  do  not  form  hydrates  for  the  same  reason. 


The  first  four  members  of  the  paraffinic  hydrocarbons 

have  molecular  dimensions  suitable  to  be  enclosed  in  cavities 

of  Structure  I  and  Structure  II.  Hydrates  of  methane,  ethane, 

and  other  saturated  and  unsaturated  hydrocarbons  have  been 

known  for  a  long  time.  Initial  hydrate  forming  conditions 

(4,  6,  7,  8) 

of  these  gases  have  been  studied  by  many  workers 

To  study  the  effect  on  hydrate  forming  conditions  of  one  gas 

in  the  presence  of  another  gas,  many  ternary  systems  have  been 

investigated.  Deaton  and  Frost studied  the  hydrate-forming 

conditions  in  the  ternary  systems,  methane -ethane -water  and 

methane -propane -water .  The  hydrate  forming  conditions  of 

binary  mixtures  of  methane  and  other  gases  were  studied  by 

(12,  26) 


Katz  and  co-workers 


A  common  trend  clearly  emerges 


out  of  the  results  obtained  on  the  above  systems. 

( 27 ) 

Van  Cleeff  and  Diepen  investigated  the  hydrate 


forming  conditions  of  nitrogen  and  oxygen  and  introduced  these 
two  permanent  gases  of  atmosphere  in  the  list  of  potential 
hydrate  formers.  As  nitrogen  is  one  of  the  constituents  of 
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the  natural  gas,  it  was  decided  to  study  the  initial  hydrate 
forming  conditions  of  a  ternary  system  N2-CH4-H2O.  In  the 
ternary  systems  previously  studied,  whenever  methane  was  present 
it  was  the  lighter  component  whereas  in  the  system  under 
investigation  methane  is  the  heavier  component. 

(33) 

Wilcox,  Carson  and  Katz  indicated  that  gas  hydrates 

behave  to  some  extent  as  solutions  of  gases  in  crystalline 

solids,  and  put  forward  the  concept  of  vapour-solid  equilibrium 

constants  which  can  be  effectively  used  to  predict  the  hydrate 

forming  conditions  of  multi -component  mixtures.  By  definition, 

an  equilibrium  constant  is  a  ratio  of  the  mole  fraction  of  a 

particular  component  in  the  gas  phase  to  that  in  the  hydrate 

phase  on  a  water-free  basis.  Consequently  the  equilibrium 

constant  for  the  lightest  component  in  a  multi -component  mixture 

will  always  be  greater  than  unity.  Equilibrium  constants  for 

methane  obtained  from  the  study  of  various  ternary  systems  were 

all  greater  than  unity,  since  methane  was  the  lightest  component. 

One  of  the  objectives  of  this  investigation  was  to  provide 

equilibrium  constants  for  methane  in  a  system  where  it  was  the 

heavier  component.  This  would  make  is  possible  to  extend  the 

(1)  ... 

curves  prepared  by  Carson  and  Katz  ,  depicting  the  equilibrium 
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constants  of  methane  as  a  function  of  temperature,  in  a  region 
where  the  value  of  equilibrium  constant  is  less  than  unity. 
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THE  SMALL  CAVITY  OF  STRUCTURE  I 


MAR  •  63 


THE  LARGE  CAVITY  OF  STRUCTURE  I 


FIGURE  1 
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THE  LARGE  CAVITY  OF  STRUCTURE  II 


FIGURE 
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THEORY 


Phase  Equilibria  in  Gas  Hydrates 

The  study  of  gas -hydrates  essentially  involves  the 
study  of  various  phases  in  equilibrium  at  different  conditions 
of  pressure  and  temperature.  Depending  on  the  critical 
temperature  of  the  hydrate  former,  and  the  working  temperature 
there  may  be  three  or  four  phases  in  equilibrium.  Below  the 
ice  point  the  phases  which  can  be  in  equilibrium  are  hydrate, 
ice,  vapour  and  liquid  phase  rich  in  hydrate  former.  The  phase 
which  can  be  in  equilibrium  above  the  ice  point  are  hydrate, 
water  rich  liquid,  vapour  and  a  liquid  phase  rich  in  hydrate 
former.  The  total  number  of  phases  present  at  any  instant  will 
depend  on  the  number  of  degrees  of  freedom  of  the  system,  which 
can  be  easily  predicted  from  Gibb's  phase  rule  stated  mathe¬ 
matically  as  follows? 

F  +  P  =  C  +  2 

where  F  =  number  of  degrees  of  freedom 
P  =  number  of  phrases  present 
C  =  number  of  components 

Thus  in  a  binary  system,  univariant  equilibria  allow 
the  presence  of  only  three  phases,  whereas  in  a  ternary  system 
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four  phases  may  be  present. 


The  following  generalization  is  thought  to  be  of  great 


help  in  a  better  understanding  of  the  phase  equilibria. 


(18) 


"The  addition  of  a  component  to  a  system  adds  a 
degree  of  freedom  to  all  features  represented  for 
the  system.  In  other  words,  areas  will  generate 
volumes,  lines  will  generate  surfaces,  and  points 
will  generate  lines.  These  features  will  have  one 
extra  degree  of  freedom  but  will  not  be  changed 
in  designation.  That  is,  a  point  representing 
vapour,  solid  and  liquid  in  equilibrium  for  a  one 
component  system  will  generate  a  line  in  space 
on  the  addition  of  a  second  component,  but  this 
line  will  represent  equilibrium  between  the  same 
phases.  Similarly  a  line  representing  vapour, 
liquid  equilibrium  will  generate  a  surface  which 
will  represent  a  liquid  saturated  with  vapour. 
Associated  with  each  surface  there  must  be  another 
surface  representing  vapour  saturated  with  a  liquid". 


Termination  of  a  phase-equilibrium  is  brought  about 
by  the  addition  of  a  phase,  by  the  loss  of  a  component  or 
by  an  attainment  of  a  critical  point.  The  effect  brought 
about  by  an  addition  of  a  phase  or  by  the  loss  of  a  component 
is  to  decrease  the  number  of  degrees  of  freedom  of  the  system 
by  one.  Thus  in  both  the  cases  the  prevailing  phase-equilibrium 
is  terminated.  The  appearance  of  a  solid  phase  brought  about 
by  cooling  a  one  component  system  consisting  of  two  phases, 
liquid  and  vapour,  down  to  its  triple  point  terminates  the 
two  phase  equilibrium,  and  the  system  becomes  invariant.  In 


* 


-  11 


a  purely  physical  process  it  is  not  possible  for  a  closed  system 
to  lose  a  component.  A  thermodynamic  system  may  however  lose 
a  component  in  a  physico-chemical  process  under  the  influence 
of  the  external  variables,  with  the  loss  of  a  degree  of 
freedom  which  eventually  leads  to  the  termination  of  the 
particular  phase  equilibrium.  Termination  of  the  phase  equili¬ 
brium  by  an  attainment  of  a  critical  point  is  very  common  in 
phase-equilibrium  studies.  Due  to  the  loss  of  identity  of 
the  two  distinct  phases  -  vapour  and  liquid  -  at  the  critical 
point  the  prevailing  vapour-liquid  equilibrium  is  terminated. 

Phase  Equilibria  in  a  Binary  System 

Figure  3  is  a  phase  equilibrium  diagram  of  a  binary 
system  consisting  of  a  hydrate  former  and  water.  Lines  BC 
and  BF  indicate  the  univariant  equilibria  with  three  phases 
in  equilibrium.  Along  BC,  three  phases  -  hydrate,  water- 
rich  layer  and  gas  -  are  in  equilibrium.  No  hydrate  will 
exist  below  BC  and  BF,  while  above  BC  and  BF,  hydrate  and 
excess  phase  will  be  present.  Figure  3  is  a  typical  pressure- 
temperature  diagram  for  a  hydrate  former  whose  critical 
temperature  is  considerably  below  the  working  temperature. 
Methane-water,  and  ethylene-water  are  typical  of  systems 
exhibiting  the  phase  behaviour  indicated  in  this  figure. 


FIG-3 
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SCHEMATIC  PHASE-DIAGRAM  OF  HYDRATES  EXHIBITING 
NO  TEMPERATURE  LIMIT 
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Figure  4  is  a  binary  phase  diagram  for  a  hydrate -former 
whose  critical  temperature  is  considerably  above  the  freezing 
point  of  water.  Lines  BC  and  BF  represent  the  same  phases  in 
equilibrium  as  in  Figure  3.  The  point  C  is  of  particular 
significance  as  it  is  invariant. 

Four  phases  in  equilibrium  at  C  are  hydrate,  water- 
rich  layer,  hydrocarbon  rich  layer  and  gas.  The  three  phases 
in  equilibrium  along  Cl  are  hydrate,  water  rich  layer,  and 
hydrocarbon  rich  layer.  Ethane -water ,  and  propane -water ,  are 
typical  systems  exhibiting  the  phase  behaviour  indicated  in 
Figure  4. 

Phase  Equilibria  in  Ternary  Systems 

The  addition  of  a  component  leads  to  an  additional 
degree  of  freedom  for  the  system.  A  ternary  system  has  one 
degree  of  freedom  more  than  a  binary  system.  This  additional 
degree  of  freedom  is  the  variable  composition  of  the  system. 
Depending  on  the  total  composition  of  the  system  the  hydrate 
will  meltover  a  range  of  temperatures  at  a  fixed  pressure. 
Similarly  there  will  be  a  range  of  decomposition  pressures 
at  a  fixed  temperature. 
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SCHEMATIC  PHASE-DIAGRAM  OF  HYDRATES  EXHIBITING 
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Quadruple  points  in  a  ternary  system  are  univariant 
and  not  invariant  as  in  the  case  of  a  binary.  There  exists 
an  unique  quadruple  point  for  every  composition.  As  a 
result  there  is  a  line  of  four  phase  equilibrium  instead  of 
a  single  point  as  in  the  case  of  a  binary. 

Hydrate  forming  conditions  in  a  ternary  system  can 
be  studied  either  by  keeping  the  composition  or  temperature 
fixed.  In  a  ternary  system,  the  three  phase  equilibrium  is 
univariant  at  constant  composition.  As  a  result,  a  family 
of  pressure-temperature  curves  each  at  a  constant  composition 
reveals  the  general  hydrate  forming  tendency  of  the  system. 

Ternary  systems  may  also  be  studied  by  determining  the  iso¬ 
thermal  pressure-composition  diagrams,  which  elucidate  the 
general  phase  behaviour  quite  clearly. 

Statistical  Mechanical  Theory  Applied  to  Hydrate  Formation 

/  OD) 

Van  der  Waals  and  Platteeuw' '  developed  the  solid 
solution  theory  for  hydrates  on  the  basis  of  statistical 
mechanical  principles.  Thermodynamic  properties  of  hydrates 
were  derived  from  a  simple  model  which  corresponds  to  the 
three  dimensional  generalization  of  ideal  localized  adsorption. 

The  structural  frame  formed  by  water  molecules  by  hydrogen  bonding 
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is  stabilized  only  by  the  gas  molecules  occupying  the 
cavities.  The  empty  host  lattice  is  metastable  and  does 
not  exist  as  such.  The  ordinary  crystal  lattice  (liquid 
water  or  ice)  is  termed  theC'-.  modification  and  the  metastable 
lattice  is  called  .  The  mathematical  analysis  presented  in 
this  chapter  is  based  on  the  following  assumptions. 

(1)  The  contribution  of  the  water  molecules  to  the  free 
energy  is  independent  of  the  mode  of  occupation  of  the  cavities. 
This  excludes  thecase  of  large  molecules  which  might  disturb 
the  lattice. 

(2)  The  mutual  interaction  of  the  solute  molecules  is 
neglected,  i.e.  the  partition  function  of  a  solute  molecule  in 
the  cage  is  independent  of  the  number  and  types  of  the  solute 
molecules  present. 

(3)  The  encaged  molecules  are  localized  in  the  cavities 
and  a  cavity  can  never  hold  more  than  one  solute  molecule. 

(4)  Classical  statistics  are  valid. 

The  behaviour  of  the  solute  molecules  in  the  cavity 
is  represented  satisfactorily  by  the  four  assumptions  given 
above.  The  objection  raised  against  the  secondassumption  is 
waived  on  the  basis  of  the  information  obtained  from  the 


numerical  calculations  and  experimental  evidence.  The  solute 
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solute  interaction  to  the  configurational  energy  is  at  the  most 
a  few  percentage  of  the  energy  of  binding  of  the  solute  mole¬ 
cules  in  the  cages. 

Partition  Function 

Mathematically  a  partition  function  is  given  by 
the  following  relation: 


P.F.  =  Q  =  g  e  i/kT  (1) 

o 

where 

gs  =  degenerancy  of  the  energy  state  i 
i  =  energy  in  the  state  i 
k  =  Boltzmann's  constant 
T  =  absolute  temperature 


Following  are  the  relations  connecting  partition  function 
with  the  thermodynamic  properties. 

F  =  -  RT  In  +  1  -  V  - — —  -  (la) 

N  V  m 


H 


whore 


_ _  In  Q  .  In  Q 

RT  - —  + - — 

In  T  In  V  ip 

F  =  free  energy/mole 
H  =  enthalpy/mole 
V  =  volume  of  the  assembly 
N  =  Avogadro's  constant 


(2) 
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Depending  on  the  relations  of  the  system  with  the 
surroundings,  it  may  be  necessary  to  use  one  of  the  three 
kinds  of  partition  functions  given  belows 

(1)  Micro  canonical 

(2)  Canonical 

(3)  Grand  canonical 


The  microcanonical  partition  function  or  some¬ 
times  known  as  the  system  partition  function  is  the  summation 
represented  in  Equation  1  for  an  individual  particle.  A 
system  in  statistical  mechanics  is  used  to  describe  microscopic 
entities  like  atoms,  molecules,  ions  and  photons. 


The  Canonical  partition  function  (sometimes  known 
as  assembly  partition  function)  is  a  multiple  of  the  partition 
functions  for  the  individual  systems. 

<P-F-> assembly  =  qN 

An  assembly  in  statistical  mechanics  consists  of  a  large  number 
of  systems,  and  it  corresponds  to  a  thermodynamic  system. 


The  grand  canonical  partition  function  (usually 
abbreviated  to  G.P.F.)  is  given  by  the  following  relations. 


^  (P.F.)  A 
(N) 


N 


(3) 


G.P.F. 


(4) 


. 
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where  , 


x 


Assembly  partition  function 
Absolute  activity  of  a  system  in 
assembly. 


Application  of  the  Theory  of  G,P0Fa  to  Clathrates 


A  clathrate  consists  of  a  cage  forming  substance  and 
a  solute.  The  partition  function  of  a  clathrate  can  be  repre¬ 
sented  as  a  multiple  of  two  part  ition  functions  as  follows  2 


(P-P‘)  clathrate  =  <P-P*>Q  X  (P'FjJ 


where 


(P.Fo)q  =  partition  function  of  a  cage  forming 

substance 


(P.F.)t  ®  partition  function  of  the  enclosed 


molecules 


(P.F.)q  =  exp  (■ 


( P . F. ) T  =  TT ( 


A 

F  /kT ) 


V  i  nq 


i  (V.nq  -  ZHTi:)  :TrNji  • 


TT  hT-  \Ji» 

J  Jl  ) 


Jl  T 

J  J 


exp 


(-  F  AT)  ^ 


NJi 


(Vj  nQ}  i 


(V.Nq  Nji)  :TrNji; 

J  J 


Nt-J 

Tt£ 1  NJi 
jnjiAJ 


(5) 


G.P.F. 


(6) 
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where 

=  number  of  cavities  of  type  i  per  molecule  of  Q 
hj  =  partition  function  of  a  J  molecule  when  encaged 
in  a  type  i  cavity. 


Equation  6,  after  summing  with  the  help  of  a  multinomial 


theorem,  takes  the  following  form 


G.P.F.  =  exp  (-  FVkT)  1T  (1  +  2  hji  ^  J) 


yiNQ 


(7) 


Composition  of  the  clathrate  is  given  by  the  following  relation 


Nk  “  ^  Nki 


A  a 

K  (  -2. 


In  G.P.F.  \  \ 

“  jT,  V,  Nq,  Aj 


b 


(8) 


Differentiating  Equation  7  with  respect  to  "^k  and 
substituting  the  result  in  Equation  8,  the  following  equation 


is  obtained? 


yjc. 


Ntfyz 


h  Ki  ^  K 


(1  + 

J 


J) 


(9) 


where 


A 


K 


kT 


^  and  ^k(T)  is  the  molecular 


K 


partition  function  of  gaseous  K  with  the  volume  factor  removed 


YKi  = 


cKi *  PK 


(1  +  *Z.  Cj^*  Pj) 


Now 


(10) 
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where 


°ki 


=  hki  ^  k 
Pk 


(ioa; 


Therefore, 


2  YKi 

K. 


=  (1  +  XCj  Pj) 


(n: 


Differentiation  of  Equation  6  with  respect  to  Nq  generates 


the  following  equations 


/  ^  In  Go  P  „  Fa  \  v 

\  ”  /  m  tr 


b  N, 


Q 


T  V 

x  0  v  0  -j 


4 

kT 


£> 


X  V •  In  ( 1  +  ^  h  T.i  \  i 

i  a  i  j . 

l  J 


where 


6 


Q 


^nQ 


is  the  chemical  potential 


of  V  modification. 


(12; 


Substituting  (1  +  2  Cj.  pj)  ^  for  (1  +  2 hj^j) 


in  Equation  12,  one  obtains 


/t 


Q 


ft 

Q 


kT 


2Y 


In  (1  +  ZCj^j) 


(13) 


Equations  11  and  13  are  very  important  relations 
connecting  vapour  pressure,  composition  and  chemical  potential 
of  the  solvent  in  a  clathrate , obeying  the  present  model. 


The  careful  inspection  of  Equation  10  reveals  the 
similarity  it  holds  with  the  Langmuir  isotherm.  The  Langmuir 


adsorption  isotherm  may  be  represented  as  follows 
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1  +  c.y  p 

where 

c-^,  are  constants,  p  is  total  pressure  of 
the  gas  phase  and  <D  is  the  fraction  of  the  surface  area  covered 
by  adsorbed  gas  molecules,. 


Equation  10,  when  written  for  a  single  component  gas 
phase,  takes  the  following  form. 


Yi 


C±  p 


1  +  C±P 


The  langmuir  isotherm  has  been  applied  to  a  number  of  cases  of 
adsorption  on  solid  surfaces,  however  it  has  its  own  limitations. 


The  constant  '  G^.  '  will  be  referred  to  as  a  L  angmuir  constant 


Evaluation  of  Langmuir  Constants 

As  shown  earlier  by  Equations  9  and  10,  the  Langmuir 
constant  has  been  defined  as  a  function  of  hj^0  Evaluation 
of  hj^  requires  a  method  to  calculate  the  potential  energy  of 
a  solute  molecule  in  the  cavity.  In  addition  t©  translational 
and  rotational  energy  of  the  molecule  there  is  an  interaction 
energy  which  should  be  included  in  the  partition  function. 
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The  theory  introduced  by  Lennar d- Jones  and  Devonshire 
for  the  study  of  liquids  provides  a  powerful  method  for  the 
quantitative  evaluation  of  the  partition  function  of  a  solute 
molecule  within  its  cavity. 

The  following  two  assumptions  in  addition  to  the  four 
assumptions  cited  earlier,  are  implied  in  the  L-J-D  theory s 

(5)  The  solute  molecules  can  rotate  freely  in  their 
cavities.  (i.e.  The  rotational  partition  function  in  the  cavity 
is  the  same  as  that  in  the  gas  phase) . 

(6)  The  potential  energy  of  a  solute  molecule  when  at  a 
distance  r  from  the  centre  of  its  cage  is  given  by  the  spheri¬ 
cally  symmetric  potential  w(r) ,  proposed  by  Lennard- Jones  and 
Devonshire. 

The  assumptions  given  above  restrict  the  applicability 
of  this  theory  to  the  monatomic  gases  like  argon,  krypton,,  and 
xenon.  Rotation  of  diatomic  molecules  in  a  cavity  is  actually 
constrained  due  to  their  non-spherical  shape.  Hindered  rotation 
of  molecules  like  oxygen  and  nitrogen  tends  to  decrease  the 
entropy  and  as  a  result  the  vapour  pressure  of  such  molecules 
is  raised. 


The  Lennard- Jones  6~12  potential  evaluates  the  average 
contribution  to  the  potential  energy  due  to  the  interaction 
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of  a  solute  molecule  with  any  of  the  elements  constituting 
the  wall  of  its  cage.  The  energy  of  interaction  is  given  by 
the  following  equations 

0(R)  =  4  £ 

(11) 

Lennard“ Jones-Devonshire v  summed  the  interactions  over  Z 

molecules  distributed  over  the  surface  of  a  sphere  and  so 
obtained  a  function  w(r)  describing  the  resulting  field „ 


(-4-) 12 

J.  R 


‘  R 


Within  the  sphere  averaged  over  all  orientations. 


w(o) 


(14) 


w(r)  -  w(o)  =  Z  ^4l(y)  ~  c<2  ->2m(y)J  (15) 

1 (y )  =  (y4  +  12y3  +  25. 2y2  +  12y  +  1)/(1  -  y) 10~1 
m(y)  =  (1  +  y)/(l  -  y) 4  -  1 

where 

y  =  (r/a)  ,  &  are  the  molecular  parameters, 

and  cK  is  a  parameter  which  measures  the  cell  radius  'a5  in 
terms  of  molecular  diameter  6  .  Mathematically 


s/~2 


With  the  assumptions  (5)  and  (6)  the  cell  partition 


25 


function  can  be  represented  as. 


=  0 (T)  exp 

r_  w(o)  |  exp  [ 

w(r)  -  w(o) 

L  RT  cel^  L 

kT  J 

1 

(16) 


0(T)  represents  the  contributions  due  to  trans¬ 
lational  and  internal  degrees  of  freedom  of  the  solute  molecule. 
The  second  term  represents  the  potential  energy  w(o)  of  the  solute 
molecule  at  the  centre  of  its  cage  referred  to  the  perfect  gasf 
and  the  integral  is  the  free  volume  of  the  solute  molecule 
wandering  in  the  cavity.  In  order  to  conform  with  the  customary 
notation  of  the  L-J-D  theory  the  free  volume  may  be  written  as 

3  i  1 

2  A  a  g,  where  the  function  g  stands  for  the  dimensionless 
integral , 

’-/•*>  -if  [-  *  -^f-]  «  <* 

Hence, 

h  =  <${T)  exp  w(o)/kTj  2*a3g 

Substituting  the  value  of  h  in  10(A) s 


2*a  • 3  gKi  exp  (- 
kT  1  ^  kT 


(18 


Equation  18  can  be  successfully  used  to  evaluate  the  Langmuir 
constants.  Numerical  integration  of  (17)  will  yield  the  value 
of  gK^  which  may  be  subsequently  used  in  (1.8)  to  obtain  CK^. 
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EXPERIMENTAL  PROGRAMME 


Equipment 


Equipment  used  in  this  investigation  was  similar  to 
that  used  by  earlier  workers  for  the  study  of  gas-hydrates  and 
volumetric  studies  at  the  University  of  Alberta.  The  equilibrium 
cell  used  in  this  work  was  the  same  as  that  used  by  Snell ^21) 
for  the  study  of  gas-hydrates. 


The  main  parts  of  the  equipment  are  as  follows; 

(1)  Equilibrium  cell 

(2)  Mercury  pump 

(3)  Rotary  seals 

(4)  Constant  temperature  bath 


(5)  Agitating  mechanism 


Equilibrium  Cell 

The  equilibrium  cell  was  a  Penberthy  liquid  level 
gauge  capable  of  withstanding  6000  psia  at  250°F,  with  a  test 
pressure  of  9000  psia.  The  cell  was  constructed  of  type  316 
stainless  steel  with  sight  glasses  on  both  sides  and  had  an 
approximate  volume  of  120  cc. 


SCHEMATIC  DIAGRAM  OF  EQUIPMENT 
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The  cell  was  enclosed  in  an  8"  O.D.  lucite  jacket 
with  lucite  discs  and  O-rings  on  both  ends  to  make  the  seal. 

Both  the  ends  of  the  cell  were  connected  with  aminco  tubing, 
which  rested  in  a  bearing  placed  in  rotary  seals. 

Two  sight  glasses,  used  on  each  side  of  the  cell  served 
as  glass  windows  for  visual  observations  of  the  phase  behaviour, 
and  were  capable  of  withstanding  9000  psia  at  100°Fo  Figure  8 
shows  the  complete  assembly  of  the  equilibrium  cell  with  the 
glasses . 

Mercury  Pump 

A  high  pressure  positive  displacement  pump,  with  a 
pressure  rating  of  8,000  psi  and  displacement  volume  of  265  cc0 
was  used  to  raise  or  lower  the  pressure  in  the  cell  by  injecting 
or  withdrawing  mercury.  Graduations  on  the  pump  indicated  the 
amount  of  mercury  injected  or  withdrawn.  Volumetric  calibration 
of  the  cell  was  performed  using  the  pump. 

Rotary  Seals 

( 21) 

Rotary  seals  similar  to  those  used  by  Snell  were 

employed  with  great  success  to  prevent  the  leakage  of  the  charge 
during  rocking  of  the  cell.  Details  of  the  seals  are  given  in 
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THE  EQUILIBRIUM  CELL 


FIGURE  6 
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Figure  7.  The  Aminco  tubes  connecting  the  equilibrium  cell 
were  made  to  pass  through  the  hollow  brass  shafts  which  rested 
in  bearings  in  the  aluminum  housings.  Synthetic  rubber  seals 
were  used  in  the  housings  to  prevent  the  leakage  of  the  coolant. 

An  Aminco  tee  with  a  V  hole  drilled  through  it  was  used  to 

make  the  seal.  The  Aminco  tube  passed  through  this  hole  and 
rotated  as  an  axle  in  bearing.  The  end  of  the  Aminco  tube  was 
welded.  A  1/16"  hole  was  drilled  through  the  side  into  the 
inside  of  the  tubing  for  the  introduction  of  the  charge  into 
the  cell. 

There  was  enough  clearance  between  the  tubing  and 
the  block  so  that  mercury  was  around  the  tubing  and  the  pressure 

therefore  transmitted  to  the  cell  at  all  times.  Teflon  0-rings 

were  tightened  down  by  steel  sleeves  enclosed  in  Aminco  nuts 
to  make  the  seal. 

Constant  Temperature  Bath 

The  contents  of  the  equilibrium  cell  were  maintained 
at  a  constant  temperature  by  circulating  varsol.  Varsol  was 
circulated  through  the  lucite  jacket  around  the  cell  using  a 
centrifugal  pump  immersed  in  an  insulated  tank  filled  with  varsol. 
A  freon  refrigeration  unit  and  two  heaters  of  1000  watt  and 
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500  watt  capacities  were  employed.  The  temperature  in  the 
bath  was  controlled  by  an  American  Instrument  Co.  bimetallic 
strip  controller.  This  controller  was  connected  through  a  relay 
and  transformer  to  the  refrigerator.  Temperature  control  was 
within  ±0.1°F. 

The  temperature  in  the  equilibrium  cell  was  measured 
with  the  help  of  an  iron-const  ant  an  thermocouple,,  inserted  through 
the  side  of  the  cell. 

Agitating  Mechanism 

A  1/6  horse-power  electric  motor  was  employed  to  drive 
an  eccentric  which  was  connected  to  the  cell  through  a  connecting 
rod.  To  provide  necessary  agitation,  the  cell  was  rocked  through 
a  small  arc  of  30°  -  35°  when  the  contents  of  the  cell  moved 
back  and  forth. 

Valves  and  Tubings 

Aminco  and  Autoclave  tubings  were  used  in  the  experi¬ 
mental  equipment.  Aminco  tubings  supplied  by  American  Instrument 
Company,  with  an  inside  diameter  of  3/32"  and  a  wall  thickness 

of  1/8"  were  used.  Both  Aminco  fittings  and  tubings  were  of 

,  o 

type  304  stainless  steel  and  rated  at  30,000  psia  at  100  F. 
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FIG.  7 
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THE  ASSEMBLY  OF  THE  EQUIPMENT 


FIGURE  8 
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Autoclave  tubings,  made  by  Autoclave  Engineers  Inc.  had  an  out¬ 
side  diameter  of  1/8"  with  a  wall  thickness  of  1/32".  The 
tubings  and  valves  were  rated  for  15,000  psia  at  100°F  and 
were  constructed  of  type  316  stainless  steel. 

The  methane  gas  used  in  this  investigation  was  obtained 
from  the  Matheson  Company,  with  a  specified  purity  of  99.8%. 

The  nitrogen  gas  was  obtained  from  the  Liquid  Carbonic 

Canadian  Corporation  Limited,  with  a  specified  purity  of  over  99.6%. 

Exp e r imental  T echn i gu e 

Figure  5  is  the  schematic  diagram  of  the  assembly 
of  equipment.  Initially  the  equilibrium  cell  was  filled  with 
mercury  with  the  help  of  the  mercury  pump.  Keeping  valves  2 
and  3  open,  some  mercury  was  collected  in  the  trap,  so  as  to  make 
sure  that  the  tubing  was  filled  with  mercury  up  to  valve  2. 

Valve  2  was  then  closed  and  the  line  from  valve  6  to  valve  1 
was  evacuated  (keeping  valve  6  closed) .  Valve  1  was  closed 
and  valve  7  was  opened  for  the  introduction  of  water  into  the 
cell.  Valves  2  and  4  were  opened,  and  by  slowly  withdrawing 
mercury  from  the  equilibrium  cell,  water  was  introduced. 
Approximately  15  -  20  cc  of  water  were  taken  each  time.  Valves 
7  and  2  were  closed  and  the  line  connecting  the  gas  cylinder 
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and  the  cell  was  evacuated.  Valve  1  was  then  closed  and  the 
desired  amount  of  gas  was  taken  into  the  equilibrium  cell.  The 
same  process  was  repeated  to  introduce  another  gas  into  the 
cell  to  make  the  desired  composition  of  the  gaseous  mixture. 

The  gas  from  a  cylinder  with  a  lower  pressure  was  always  taken 
in  first.  Any  predetermined  mixture  of  the  gases  could  be  pre¬ 
pared  by  adjusting  the  initial  pressure  of  the  pure  gas  samples. 
The  charge  thus  taken  into  the  cell  was  allowed  to  mix  completely 
by  rocking  the  cell  for  15  -  20  minutes.  Special  care  was  taken 
to  see  that  valves  2  and  4  were  kept  closed  during  the  experi¬ 
ment.  Later  a  small  sample  of  this  charge  was  chromatographically 
analysed . 

Technique  for  the  Determination  of  Initial  Hydrate- Forming 

Conditions 

Once  the  composition  is  fixed ,  a  ternary  system 
behaves  exactly  like  a  binary  system  as  regards  its  initial 
hydrate- forming  conditions.  Therefore,  the  procedure  followed 
for  the  determination  of  initial  hydrate  forming  conditions 
for  a  ternary  system  with  a  fixed  composition  is  exactly  the 
same  as  that  for  a  binary  system.  At  a  particular  temperature 
the  equilibrium  pressure  of  the  three  phases  can  be  approximately 
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determined.  The  contents  of  the  cell  were  kept  at  a  pressure 
slightly  above  this  predetermined  pressure  and  the  system  was 
subjected  to  supercooling  of  8°F  to  12°F.  The  cell  was 
allowed  to  rock  for  30  to  60  minutes,  when  the  hydrate  phase 
usually  appeared.  Formation  of  hydrate  was  accompanied 
with  a  drop  in  pressure.  The  temperature  of  the  system  was 
slowly  raised  to  decompose  all  but  a  few  crystals  of  hydrates. 
With  the  rise  in  temperature,  the  pressure  also  rose  slowly 
due  to  decomposition  of  the  hydrate.  When  very  few  crystals 
of  the  hydrate  were  seen  adhering  to  the  glass  wall,  the 
temperature  and  pressure  were  noted.  The  final  pressure  of 
the  system  was  very  close  to  the  initial  pressure.  At  this 
stage,  lowering  the  temperature  by  half  a  degree  to  one  degree 
was  enough  to  obtain  sufficient  hydrate  crystals.  The  pressure 
of  the  system  was  raised  and  the  same  procedure  was  repeated 
to  get  another  point. 

Technique  for  the  Determination  of  Pressure  -  Composition 

Diagram  for  N? --CH4-H ?Q  System 

In  a  ternary  system,  the  composition  of  the  two 
phases,  gas  and  hydrate,  in  equilibrium  with  the  water  phase, 
will  depend  on  the  pressure  and  temperature  of  the  system.  The 


. 
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method  followed  to  obtain  data  on  this  aspect  of  the  three  phase 
equilibrium  is  briefly  described  below. 

Initially  the  charge  consisting  of  water  and  desired 
composition  of  the  gas  was  taken  in  the  cell  as  described  earlier. 
The  charge  was  compressed  to  a  pressure  which  was  25%  -  30% 
greater  than  that  needed  for  initial  hydrate  forming  conditions. 
The  temperature  was  set  at  the  desired  value  and  after  closing 
valve  2,  the  cell  was  allowed  to  rock  for  four  to  six  hours. 
Special  care  was  taken  to  see  that  the  plunger  of  the  mercury 
pump  was  maintained  in  its  position  throughout  the  experiment. 

In  the  end  when  equilibrium  between  the  phases  was  established 
the  pressure  was  constant.  The  drop  in  pressure  was  a  direct 
function  of  the  amount  of  hydrate  crystals  formed.  After 
noting  the  final  pressure  the  sample  of  the  gaseous  phase  was 
analysed  chromatographically . 

Hydrate  composition  was  determined  by  direct  analysis 
as  well  as  by  calculations.  First  the  gas  from  the  equilibrium 
cell  was  isobarically  transferred  to  the  auxiliary  cell.  Hydrate 
crystals  were  retained  at  the  top.  Drops  of  water  and  mercury 
falling  into  the  auxiliary  cell  indicated  the  complete  trans¬ 
fer  of  gas  from  the  equilibrium  cell  to  the  auxiliary  cell. 

Valve  2  was  closed  and  the  lines  connecting  the  equilibrium  cell 
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to  the  auxiliary  cell  were  evacuated  to  prevent  contamination 
of  the  sample  with  the  traces  of  the  gas  remaining  in  the  line. 
The  pressure  in  the  cell  was  decreased  and  the  temperature 
was  raised  to  decompose  the  hydrate.  Later,  the  sample  of 
the  gas  generated  on  decomposition  of  the  hydrate  was  chromato- 
graphically  analysed.  An  indirect  method  of  determining  the 
hydrate  composition  from  the  gas-phase  compositions  and  pressure 
data  is  given  elsewhere. 

The  procedure  described  above  was  repeated  to  cover 
the  range  of  pressures  in  which  the  two  gases  form  hydrates. 
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EXPERIMENTAL  RESULTS 


Initial  Hydrate-Forming  Conditions 


Six  mixtures  of  nitrogen  and  methane  were  studied 
to  obtain  data  on  equilibrium  pressures  at  various  temperatures 
for  the  three  phase-HLG-equilibrium.  The  temperature  and 

pressure  ranged  from  32°F  to  72°F/  and  300  psia  to  5,100  psia 
respectively. 

The  pressure  temperature  data  are  presented  in  Table 
1.  The  pressure  temperature  data  for  hydrates  of  pure  nitrogen 
and  methane  were  also  obtained  and  found  to  agree  very  well 
with  the  data  of  Van  Cleeff  and  Diepen^^,  otto^^),  and 
Snell ( 21) . 


Figure  9  represents  the  hydrate  forming  conditions 
of  pure  nitrogen,  pure  methane  and  mixtures  of  nitrogen  and  meth¬ 
ane.  The  pressure  temperature  loci  in  this  figure  show  a  general 
trend  which  one  might  expect  from  the  published  data  on  other 
ternary  systems.  Close  observation  of  the  family  of  plots 
indicates  the  interesting  fact  that  the  slope  of  the  pressure- 
temperature  curves  increases  with  the  rise  in  temperature  and 
it  increases  at  a  faster  rate  for  a  gaseous  mixture  with  lower 
nitrogen  content.  At  a  certain  equilibrium  pressure,  the 
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FIG.  10 

ISOTHERMAL  DATA  ON  INITIAL  HYDRATE  FORMING 
CONDITIONS  FOR  N2-CH4-H20  SYSTEM 
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ISOBARIC  DATA  ON  INITIAL  HYDRATE  FORMING 
CONDITIONS  OF  N2-CH4-H20  SYSTEM 
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rate  of  change  of  equilibrium  temperature  with  the  percentage 
of  nitrogen  in  gas  containing  more  than  50%  N2  is  higher  than 
that  existing  for  gases  containing  less  than  50%  N2.  Conse¬ 
quently  pressure  temperature  curves  for  gaseous  mixtures  with 
a  fixed  difference  in  nitrogen  content  are  closely  spaced  as 
the  nitrogen  content  decreases. 

It  is  interesting  to  note  that  the  equilibrium  pressures 
for  the  mixture  with  87 . 1%  CH^  are  lower  than  those  for  pure 
methane  above  67°F.  A  similar  phenomenon  was  observed  by 
Otto^^)  and  Snell  (21)  during  the  investigation  of  CH^-^H^-H^O 
system.  Equilibrium  pressures  for  mixtures  of  methane  and 
ethylene  were  lower  than  those  for  pure  ethylene  above  a  cer¬ 
tain  temperature. 

Pressure  Composition  Diagrams 


The  pressure  composition  diagrams  for  N2-CH4-H2O 
systems  at  32°F,  39.7°F,  and  44°F  were  determined  experimen¬ 
tally.  They  are  shown  in  Figures  12,  13,  and  14A. 

It  is  apparent  from  the  diagrams  that  the  points 
indicating  the  gas-phase  compositions  in  the  diagrams  do  not 
show  as  much  deviation  from  the  mean  curve  as  the  points  indica 
ting  the  hydrate  phase  compositions.  This  discrepancy  for  the 
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FIG.  12 


PRESSURE  COMPOSITION  DIAGRAM  FOR 
N2-CH4-H20  SYSTEM  AT  32  °F 
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hydrate  branch  may  be  attributed  to  the  fact  that  during  the 
transfer  of  gas  from  the  equilibrium  cell  to  the  auxiliary  cell, 
truly  isobar ic  conditions  were  not  maintained  and  the  hydrate 
composition  might  have  changed  due  to  disturbed  equilibrium 
conditions.  The  likelihood  of  occlusion  of  gas  in  the  hydrate 
crystals  appears  to  be  another  plausible  explanation  for  this 
discrepancy.  The  calculated  values  of  hydrate  compositions  do 
not  show  very  much  deviation  from  the  experimental  values.  In 
the  worst  cases  the  deviation  is  less  than  15%.  It  is  inter¬ 
esting  to  note  that  in  most  of  the  cases  the  calculated 
values  of  hydrate  compositions  give  higher  methane  concen¬ 
tration.  As  the  indirect  method  for  the  calculation  of  hydrate 
composition  as  described  in  Appendix  B  is  very  sensitive  to  the 
composition  and  compressibility  factors  of  the  gaseous  mixtures, 
the  experimental  values  of  the  hydrate  compositions  were  thought 
to  be  more  reliable  and  have  been  used  in  the  above  diagrams. 
Compressibility  factors  used  in  these  calculations  were  com¬ 
puted  from  an  equation  of  state  and  not  obtained  experimentally. 
The  indirect  method  for  the  determination  of  the  hydrate  phase 
composition  is  illustrated  in  Appendix  B. 
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THEORETICAL  CAL  CUT..  AT  IONS 


Prediction  of  Pressure  -  Composition  Diagrams 


One  of  the  objectives  of  this  investigation  was  to 
apply  the  solid  solution  theory  to  the  experimental  results. 
It  was  decided  to  obtain  the  pressure-composition  diagrams 
at  32°F  for  three  ternary  systems,  N2-CH4-H2O,  CH4-C2H5-H2O 
and  CH4-C2H4-H2O ♦ 

The  underlying  assumptions  and  the  various  steps 
involved  in  these  calculations  are  given  below. 


At  equilibrium,  three  phases  -  hydrate,  water  and  gas 
coexist.  For  a  hydrate  to  be  stable  with  respect  to 
modification  of  v/ater,  (ice  or  liquid  water)  ; 


4  1 

1  Q  Q 


(19) 


writing 


At 


(T,P) 


'  Q  "  Q 


cx. 


(20) 


for  the  difference  between  the  chemical  potentials  of  the  two 
modifications,  which  is  a  function  of  temperature  and  total 
pressure,  and  substitution  oj.  Equation  13  into  Equation  30 


above,  gives, 
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2  4  In  (1  +  T.  cKi  PK)  =  AH-/RT 


I< 


(21) 


For  an  equilibrium  hydrate  the  sign  of  equality  in  the 
Equation  21  holds .  Nitrogen  and  methane,  being  small  molecules 
relative  to  CgHg  and  SFg  occupy  both  the  cavities  of  Structure  I 
and  the  Equation  21  will  have  the  following  form  for  the 
Ng-CH^-HgO  system. 

^  In  (1  +  Cm  •  x  •  p  +  CN  .  — -  )  +  '  ?  In  (1  +  CM 

-*-1  1  l-x»p  z 


x  *  p  +  CN 


2  l-xep 


_A_T 

)  “  RT 


(22) 


The  values  of  the  various  constants  needed  in  the 
above  equation  are  as  follows^29); 


Y  _ 


1/23 


=  0.167  heal. /mol 


y2  =  3/23 


Langmuir  Constants 


1 


Nitrogen 

0.0474 

0 . 0606 


Methane  Ethane 

0.2361  4.67 


0.2685  4.350 


Ethylene 

2.185 

14.57 


Langmuir  constants  were  evaluated  from  the  data  on  fractional 

occupation  and  equilibrium  pressures  reported  by  Van  der  Waals 

(30) 


and  Platteeuw 
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Writing 


cMl  #  X  +  CNi  #  (1  -  x)  =  D± 


and  CM^  •  x  +  •  ( 1  -  x)  =  D 


N- 


Equation  22  takes  the  following  form, 


(1  +  DXP)  (1  +  D2P)3  =  e23(A^/RT) 


On  expanding, 


2  r>2 


1  +  (3  D2P  +  DXP)  +  3  D2^P^  +  3  DxD2  P“ 


Writing 

(a) 


(b) 


+  D23P3  +  3  DxD22P3  +  DiD23  p4  =  e23 


3  D2  +  Di 


D1  *  D2 


=  4  C 


D23  =  3  DXD22 


D1  *  D2 


=  4  a 


(c) 


3  P2^  +  3  DiD2  =  6  b 

D1  •  D23 


and  (d) 


1  _  e23(A/y RT)  _ 
D1  *  d23 


Substituting  these  values  in  Equation  23,  the 
following  quartic  equation  is  obtained 


p4  +  4  ap  +  6  bp  +  4  pc  +  d  -  0 


(23) 


(  A/4  /RT) 


(24) 
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Roots  of  the  above  equation 

=  -a  +  */~u  + 

R2  =  -a  +Npu" 

R3  =  -a  -  sj~u*  + 


are  , 

\fv  +  ^  W 

\[  v  +  ^fw7 


\[  V  -  \Tw~ 


N 


i  V  -  fw" 


where 

u  =  g  +  1,  v  =  2g  -  1,  w  =  4u2  +  3k  -  12gl 

g  =  a^  -  b,  h  =  b^  +  c^  -  2abc  +  dg,  k  =  4/3  ac  -  b^  -  d/3 
1  =  1/2  (h  +  >|h2  +  k3)1/3  -l-  1/2  (h  -  \fh2  +  k3)  i/3 


Evaluation  of  the  roots  Ri,  R2,  R3  and  R4  was  undertaken  by  writing 
a  programme  for  the  computer  IBM  1620. 


The  method  illustrated  above  gives  the  direct  solution 
of  Equation  24  without  going  into  trial  and  error  procedures « 

Out  of  the  four  roots  only  one  of  them  gives  the 
desired  value  of  the  equilibrium  pressure  for  a  certain  gas  phase 
composition. 


Composition  of  the  hydrate  phase  was  computed  from 


the  following  equations; 

z:  Vi 


YKi 


X 


x  V  i 
I<,  i 


(25) 
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FIG.  16 

PRESSURE  COMPOSITION  DIAGRAM  FOR 
CH4-C2H4-H2O  SYSTEM  AT  32  °F 
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where 

cKi  *  Pk 

yK.  _  (26) 

1  +  2=-  i  i  •  YKi 

K 

Theoretically  predicted  pressure  composition  diagrams 
for  the  three  ternary  systems,  N2-CH4-H2O,  CH4-C2H6-H2O  and 
CH4-C2H4-H2O  are  given  in  Figure  14,  15,  and  16  respectively. 
Theoretically  predicted  diagrams  have  the  corrected  values 
of  pressures  as  an  ordinate.  Experimental  values  of  pressures 
are  higher  for  all  the  gases.  This  is  particularly  true  for 
N2,  C2H4  and  C2II5.  Theoretically  predicted  values  of  equilibrium 
pressures  were  corrected  as  follows, 

P  ,  ,  =  P,,  ,  .  ,  (r,  •  x  +  r0  (1  -  x)J 

corrected  theoretical  \  1  ^  ' 

=  (Pa)  experimental 

(PA)  theoretical 

(Pb)  experimental 

r  2  =  ‘ 

(PB)  theoretical 

and  x  ~  Mole  fraction  of  A 

rl'  r2  were  evaluated  at  the  temperature  at  which 
pressure  composition  diagrams  were  investigated. 
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FIG.  17 

EQUILIBRIUM  CONSTANTS  FOR  METHANE 
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EQUILIBRIUM  CONSTANTS  FOR  N2 
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Evaluation  of  the  Equilibrium  Constants 


As  stated  earlier,  Wilcox,  Carson  and  Katz  were  the 
first  to  propose  a  solid-solution  theory  and  develop  the  concept 
of  solid— vapour  equilibrium  constants  for  gas  hydrates. 


The  equilibrium  constant  has  been  defined  as. 


Vs  =  y/x  = 


Mole  fraction  of  gas  in  vapour  phase 
Mole  fraction  of  gas  in  hydrate  phase 


V-s  ra^ios  f°r  niany  hydrocarbon  and  non-hydrocarbon  gases  were 
obtained  experimentally  from  the  study  of  ternary  systems  by 
earlier  workers^-'  ^2'  2^)  . 


The  pressure  composition  diagrams  obtained  by  the 
statistical  mechanical  treatment  of  the  solid-solution  theory 
to  the  gas  hydrates  also  give  the  values  of  equilibrium  constants 
at  32°F  for  N2,  CH4,  C2H g,  and  C2H4.  Figure  19  represents  the 
values  of  equilibrium  constants  for  CH4  obtained  from  the 
theoretically  predicted  pressure-composition  diagrams  for  the 
three  ternary  systems.  Experimental  values  of  the  equilibrium 
constants  for  methane  are  also  shown  in  Figure  19  for  comparison. 
The  values  of  the  equilibrium  pressures  plotted  in  this  figure 
are  the  corrected  values  of  the  pressures  theoretically  obtained. 
The  values  of  the  equilibrium  constant  for  CH4  at  pressures  of 
100  to  300  psia  were  obtained  from  the  data  of  Carson  and  Katz v 1  . 
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FIG.  19 

EQUILIBRIUM  CONSTANTS  FOR  METHANE 
AT  32  °F 
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These  values  of  equilibrium  constants  were  slightly  lower  than 
those  obtained  from  the  theory.  The  agreement  existing  between 
the  experimental  and  theoretical  values  of  equilibrium  constants, 
proves  the  ability  of  the  statistical  mechanical  treatment 
of  the  solid-solution  theory  to  predict  the  distribution  of 
a  component  among  two  phases,  gas  and  hydrate,  with  a  fair 
degree  of  accuracy.  The  maximum  deviation  is  less  than  25%. 

The  values  of  equilibrium  constants  obtained  from 
the  experimental  pressure  composition  data  at  32°F,  39. 7°F  and  44°F 
were  used  in  extending  the  equilibrium  constant  curves  for  methane 
in  a  region  where  the  value  of  the  equilibrium  constant  is  less 


than  unity. 


62 


DISCUSSION  OF  RESULTS 

Initial  Hydrate-Forming  Conditions 

Results  on  initial  hydrate- forming  conditions  for 
pure  nitrogen,,  pure  methane,  and  mixtures  of  nitrogen  and  methane 
are  depicted  in  Figure  9.  It  becomes  apparent  that  the  behavior 
of  the  ternary  system  under  investigation,  N^CH^-H-^O,  is  very 
similar  to  that  of  CH4-C2H6-H20  or  CH4-C3H8-H20  as  far  as 
initial  hydrate-forming  conditions  are  concerned.  The  same 
general  trend  is  exhibited  by  as  regards  the  effect 

brought  about  by  the  small  concentrations  of  the  heavier 
component  on  the  values  of  equilibrium  pressures.  Small  amounts 
of  methane  bring  about  a  large  drop  in  the  equilibrium  pressures 
at  the  same  temperatures.  Equilibrium  pressure  responds  to 
temperature  in  a  manner  which  is  similar  to  that  exhibited  by 
other  ternary  systems.  It  is,  however,  interesting  to  note 
that  the  slope  of  the  P-T  plot  shows  a  systematic  variation 
as  the  methane  content  of  the  gas  increases  from  zero  to  one 
hundred  percent.  The  slope  of  the  pressure-temperature  curve 
for  pure  nitrogen  is  constant  throughout  the  range  of  temperature 
studied.  For  a  particular  mixture  of  nitrogen  and  methane 
the  slope  of  the  pressure-temperature  curve  is  low  at  lower 
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temperatures,  but  rises  considerably  with  the  rise  in  tempera¬ 
tures.  Pressure-temperature  curve  for  methane  clearly  exhibits 
this  fact.  It  is  interesting  to  note  that  the  equilibrium 
pressures  for  the  mixture  with  87.1%  CH^  are  lower  than  those 
for  pure  methane  above  67°F. 


The  Clapeyron  equation  defines  the  slope  of  a  pressure 
temperature  curve  of  a  univariant  phase-equilibrium  by  the 
following  equation. 

dP  _  AH  (27) 

dT  T  AV 

The  physical  process  may  be  represented  by  the  following  equation 

G  '  n  H2°(s) - *  G(g)  +  n  H2° ( 1)  (28) 

Av  can  be  approximately  equated  to  — where  Z  is  the  compressi¬ 
bility  factor  of  the  gas.  Equation  1,  after  substituting  for 
AV,  takes  the  following  form. 

1  dP  AH 

p  dT  ZPT2 


d  In  P 
dT 


AH 

ZRT2 


(29) 


The  slope  of  a  pressure  temperature  curve  on  a  semi- 
logarithmic  plot  is  given  by  the  right  hand  side  of  Equation  3. 
For  a  small  temperature  difference  a  change  in  the  quantity  /t^ 
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may  be  assumed  to  be  very  small  (T  being  the  absolute  tempera¬ 
ture*  .  Consequently  the  slope  will  change  according  to  the 

•  »  .  Atj 

variation  in  the  quantity  ~~~  . 

Ah  at  temperatures  T  and  T  +  At  may  be  correlated 
according  to  the  following  relation  assuming  the  applicability 
of  Kirchoff“s  law* 

h(T  +  AT)  =  Ah(T)  +  (ACp)AT 

Acp  is  the  difference  in  specific  heats  of  the  products 
and  reactants.  According  to  Equation  2,  ACp  may  be  represented 
by  the  following  relation. 

ACp  =  n  •  Cp(H20)  +  Cp(G)  -  Cp  (Hydrate) 

As  Gas-Hydrate  is  a  solid  phase ,  it  can  be  intuitively  assumed 
that  Ac  is  a  positive  quantity.  Hence  it  can  be  safely  assumed 

Jr 

that  AH  rises  with  the  rise  in  temperature.  The  variation  of 
compressibility  factors  for  nitrogen  and  methane  along  the 
pressure-temperature  curves  is  shown  in  Figure  20.  The  compress i 
bility  factor  for  nitrogen  rises  with  the  rise  in  temperature 
and  that  for  methane  drops  with  the  rise  in  temperature.  Con¬ 
sequently  it  may  be  expected  that  the  overall  variation  of  the 
quantity  AH/Z  for  nitrogen  is  insignificant,  and  as  a  result  the 
slope  of  the  equilibrium  curve  is  constant.  On  the  other  hand. 


COMPRESSIBILITY  FACTOR 
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FIG.  20 
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variations  in  H  and  Z  tend  to  increase  the  value  of  the  slope 
of  the  pressure-temperature  curve  for  methane. 

The  above  discussion  can  be  successfully  applied  to 
explain  the  overlapping  of  the  pressure -temperature  curves 
for  87.1%  CH4  gas  and  pure  methane.  It  is  possible  that  at 
higher  temperatures  the  variation  in  ah/Z  for  87.1%  CH  gas 
is  less  than  that  for  the  pure  gas. 

This  also  explains  the  paradoxical  behaviour  of  the 
variable  R  in  Figure  22 ,  which  represents  the  ratio  of  the 
equilibrium  pressures  for  a  mixture  and  pure  methane,  at  the 
same  temperature.  It  was  previously  expected  that  R  would 
increase  with  the  rise  in  temperature  for  a  mixture  containing 
certain  amounts  of  nitrogen.  As  a  matter  of  fact  R  decreases 
with  the  rise  in  temperature.  The  smaller  values  of  R  for 
higher  temperatures  can  be  explained  on  the  basis  of  the  varia¬ 
tions  of  the  slope  of  pressure  temperature  curves  with  tempera¬ 
ture  for  pure  methane  and  a  gaseous  mixture. 

Pres sure -Compos it ion  Diagrams 

Statistical  mechanical  treatment  of  the  solid-solution 
theory  of  the  gas -hydrates ,  though  inadequate,  describes  the 
actual  behaviour  in  a  fairly  satisfactory  manner.  Figure  12 
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FIG.  22 

EFFECT  OF  NITROGEN  ON  THE  EQUILIBRIUM 
PRESSURE 
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compares  the  experimental  pressure-composition  data  with 
those  obtained  from  the  theory  at  32°F.  The  actual  values 
of  the  equilibrium  pressure  predicted  from  the  theory  were 
corrected  for  the  inherent  defect  of  the  potential  model  by  an 
empirical  formula.  Considering  the  limitations  of  the 
experimental  technique,  the  theory  appears  to  handle  the  actual 
situation  very  well.  The  deficiency  of  the  theory  in  pre¬ 
dicting  the  pressures  accurately  may  be  attributed  to  the 
assumption  of  spherical  shape  of  the  molecules  and  central  force 
field .  Large  non-spher  ical  molecules  of  nitrogen  and  ethane  can 
not  rotate  freely  in  the  cavities  and  it  is  thought  that  this 
hindered  rotation  is  the  cause  of  the  large  deviations  in  the 
equilibrium  pressures.  The  general  formulas  for  obtaining  the 
distribution  of  a  particular  component  in  the  gas  and  hydrate 
phases  still  remain  applicable.,  The  correction  applied  to  the 
theoretically  predicted  values  of  equilibrium  pressures  using 
an  empirical  formula  is  essentially  the  same  as  determining  the 
Langmuir  Constants  for  the  gases  from  experimental  values  of 
equilibrium  pressures. 

To  extend  the  applicability  of  the  statistical 
mechanical  theory  to  large  non-spherical  molecules,  it  is  nec¬ 
essary  to  obtain  the  Langmuir  constants  for  the  gases  from  the 
experimental  values  of  the  equilibrium  pressures.  Large 
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deviations  may  be  anticipated  if  theoretically  obtained  values 
of  Langmuir  constants  are  employed  in  the  calculations  of 
equilibrium  pressures.  At  this  stage,  the  procedure  employed 
by  Van  der  Waals  and  Platteeuw^^  to  obtain  the  modified 
values  of  Langmuir  constants  for  propane  and  methane  in  the 
ternary  system  CH4-C3H0  H2O  may  be  worthy  of  a  critical  review. 
The  mixed  hydrates  of  propane  and  methane  grow  in  Structure  II 
due  to  the  large  size  of  the  propane  molecule.  The  fractional 
occupation  of  the  cavities  (larger  cavities  of  Structure  II) 
by  SFg  and  C3H0  was  assumed  to  be  the  same.  The  fractional 
occupation  of  the  cavities  for  SF^  was  calculated  from  the 
knowledge  of  the  equilibrium  pressure  of  SF^  hydrate  at  -3°C, 
and  theoretical  value  of  the  Langmuir  constant  which  is  found 
to  be  570  atm From  this  information  it  was  possible  to  obtain 
the  value  of  ^ for  Structure  II. 

The  assumption  of  equal  fractional  occupation  of  the 
cavities  by  C3H0  and  SF6,  along  with  the  knowledge  of  Ah'  for 
Structure  II  made  it  possible  to  obtain  the  value  of  Langmuir 
constant  for  propane. 

Evaluation  of  Langmuir  constants  for  CH4  entailed 
the  adjustment  of  the  molecular  parameter  feA,  in  order  to  match 
the  experimental  value  of  the  equilibrium  pressure  at  -3°C 
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with  that  obtained  from  the  theory  using  the  new  value  of 
the  molecular  parameter,,  Langmuir  constants  for  methane  in 
Structure  II  were  obtained  using  the  new  value  of  the 
parameter  along  with  other  constants  reported  in  the  liter a- 

ture . 


As  Lennard- Jones  potential  can  handle  only  spherical 
molecules  satisfactorily,,  it  is  necessary  to  search  for  a  poten¬ 
tial  model  capable  of  handling  non-spherical  molecules ,  by  taking 
into  account  the  forces  of  attraction  correctly  without  the 
restrictive  assumption  of  the  central  force-field .  Corner's 
potential  model ^  and  Kihara's  core  models appear  quite 
promising.  Once  the  mathematical  technique  is  established  to 
handle  these  potential  models,  statistical  mechanical  treatment 
of  non-spherical  molecules  can  satisfactorily  account  for  the 
phase  behaviour  of  hydrates. 

Equilibrium  Constants 

As  stated  earlier,  statistical  mechanical  treatment 
of  gas-hydrates  using  the  Lennard-Jones  potential  satisfactorily 
predicts  the  distribution  of  a  component  in  the  gas  and  hydrate 
phases.  Equilibrium  constants  therefore  can  be  predicted 
fairly  accurately  from  the  theory.  Theoretically  predicted 


r 
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values  of  equilibrium  constants  for  methane  at  32°F  are  coin- 
pared  with  those  obtained  from  the  experimental  results  in 
Figure  19.  Agreement  among  the  two  sets  of  values  confirms 
the  ability  of  the  theory  to  predict  equilibrium  constants  with 
a  fair  degree  of  accuracy. 


OUTLINE  FOR  FUTURE  WORK 


The  application  of  the  solid  solution  theory  to  the 
phase  behaviour  of  gas  hydrates  might  include  the  following. 

(1)  A  further  test  of  the  Lennard- Jones-Devonshire  theory  in 
predicting  the  equilibrium  constants  in  a  ternary  system 
consisting  of  gases  of  relatively  small  molecules  is 
suggested.  It  might  be  necessary  to  apply  an  empirical 
correction  to  the  predicted  values  of  equilibrium  pressures. 

(2)  As  the  mixed  hydrates  of  the  ternary  system 

grow  in  Structure  I,  it  was  not  necessary  to  adjust  the 
molecular  parameters  of  nitrogen.  It  is  suggested  that  a 
ternary  system  consisting  of  nitrogen  and  a  hydrocarbon 
which  forms  hydrates  in  Structure  II  may  be  studied.  The 
theoretical  investigation  of  such  a  system  will  necessitate 
the  determination  of  the  Langmuir  constants  for  nitrogen 
in  Structure  II.  It  is  suggested  that  the  same  technique 
may  be  used  to  obtain  the  Langmuir  constants  for  relatively 
small  molecules  in  Structure  II.  The  prediction  of  the 
behaviour  of  a  multicomponent  system  is  feasible,  once 
the  Langmuir  constants  are  available. 
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(3)  The  evaluation  of  the  constant  !g^'  for  various  gases 
may  be  undertaken  with  a  view  to  compile  a  list  of  the 
Langmuir  constants  which  can  be  used  to  predict  equilibrium 
constants . 

(4)  For  a  satisfactory  explanation  of  the  behavior  of  non- 
spherical  molecules  in  cavities  it  is  necessary  to  apply 
potential  models  which  correctly  account  for  the  forces 

of  attraction  without  the  assumption  of  the  central  force- 
field.  Kihara’s  core  models  and  Corner's  potential  model 
may  be  successfully  applied  for  the  prediction  of  phase- 


behavior  of  hydrates. 
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NOMENCLATURE 


9S 

G.P.F. 
hJi'  hKi 


H 

H 


I 


k 


k 


v-s 


L 


1 


L 


2 


P 


PJ'  Pk 
Q 


Cell  radius 
Number  of  components 
Heat  Capacity 

Langmuir  constant  of  a  k  molecuie  xn  a  type  x  cavity 
Number  of  degree  of  freedom 
Free  energy/mole 

A  measure  of  the  free  volume  of  the  solute  molecule 
wandering  in  the  cavity 
Degeneracy  of  the  energy  state  i 
Grand  partition  function 

Partition  function  of  a  J  (or  K)  molecule  in  a  type  i 
cavity 

Hydrate  phase 
Enthalpy/mole 
Ice  phase 

Boltzmann Js  constant 

Equilibrium  constant 

Liquid  phase  rich  in  water 

Liquid  phase  rich  in  hydrate  former 

Number  of  phases 

Partial  pressure  of  J  or  K  component  gas 
Partition  function 
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R 

T 

V 

x 

y 

z 

z 

o< 

C< 

y. 

£. 

i 

<$K 

A 

A* 


A  gas  capable  of  hydrate  formation 

Absolute  temperature 

Volume 

Mole  fraction  of  a  gas 

Mole  fraction  of  a  gas  in  the  gas  phase 
Probability  of  finding  a  K  molecule  in  a  cavity  of 
type  i . 

Number  of  oxygen  atoms  surrounding  the  two  types  of 
cavities 

Compressibility  factor 

Stable  modification  of  water  (ice  or  liquid  water) 

A  dimensionless  parameter  of  the  cell  geometry 
Metastable  modification  of  water 

Number  of  cavities  of  the  ith  type  per  water  molecule 
Energy  of  the  state  i 

Partition  function  of  a  K  molecule  with  the  volume 
factor  removed 
Absolute  activity 
Chemical  potential 


A  molecular  parameter 
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APPENDIX  A 

Calculation  of  Compressibility  Factors  for 
Nitrogen-Methane  Mixtures 


As  no  adequate  data  on  the  compressibility  factors 
of  N2-CH4  mixtures  were  found  in  literature,,  it  was  decided 

/  g  \ 

to  use  an  equation  of  state  proposed  by  Keys  and  Burks'  •  to 
obtain  the  compressibility  factors  for  mixtures.  The  constants 
for  the  mixtures  were  obtained  by  a  linear  combination  of  the 
constants  of  pure  components. 


The  equation  of  state,  when  applied  to  mixtures  has 
the  following  form, 


P  = 


IRi  Wi 


-  7*  JiVfi 

V  -  IB 

1  V 


2  Aj^  wx  ) 

(V  +  £  liw±  )2 


(I) 


where , 

P  =  Pressure  in  atmosphere 

=  Universal  gas  constant  of  the 
ith  component  in  the  (  .&tm,  cc j 
mixture  9m*  0K 


w.  -  weight  fraction  of  the  ith  component  in 
the  mixture 

On 

T  =  Temperature  k 
V  =  Volume,  cc/gm. 


B.  <X 


constants  of  the  ith  component  in  an  equation 
of  state. 


-  80 


The  constants  for  nitrogen  and  methane  are  as  follows 


Methane 

Nitrogen 

A 

10151.40 

1650.50 

B 

3.961 

1.65 

C< 

2.857 

0.992 

1 

0.536 

0.313 

R 

5.1173 

2.9286 

constants 

for  mixtures  were 

obtained  by  using  the 

following  rules. 

^mix.  =  (  Anitrogen  x  WN2  +  ACH4  x  WCH^  } 

Bmix.  =  BN2  x  wN2  +  bCH4  x  wCH4 

X  WN2  +  ^CH4  x  wCH4 


1  =  1n2  x  wn2  +  1ch4  x  wCH4 


R 


%2  x  wN2  +  RCH4  x  wCH4 


writing. 


PRjL  and  pR2  for 


2  Rj  Wi 


and 


(z  /Aj  Wj  ) 


V  -  ZRi  e 


T  W  -i  O 

■=— i (V  +  Zl-W^  )Z 

V 


respectively,  Equation  I  takes  the  following  form: 


p  =  pRi  T  -  Pr2 


(II) 
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pv  _  PRi;v  pr2*v 

RT  R  ~RT 


(III) 


The  values  of  PR^  and  Pj^  were  computed  on  the  Computer  IBM 
1620,  covering  a  range  of  volume  for  5  to  60  cc/gm*  Using  these 


values  of  Pr^  and  Pr^*  compressibility  factors  at  any  tempera¬ 
ture  can  be  obtained. 


Vol , 

8% 

Pri 

N2 

Pr2 

20% 

PR1 

N2 

rr2 

32% 

Pri 

n2 

Pr 

R2 

5 

1.763 

302.310 

1.629 

261.072 

1.497 

222.501 

15 

0.417 

38.225 

0.388 

32.805 

0.361 

27.783 

25 

0.229 

14.136 

0.214 

12.115 

0.200 

10.246 

35 

0.157 

7.296 

0.1474 

6.250 

0.138 

5.283 

45 

0.119 

4.443 

0.112 

3.804 

0.105 

3.214 

55 

0.0961 

2.986 

9.059 

2.556 

0.085 

2.159 

44% 

N2 

56% 

N2 

68% 

N2 

5 

1.366  186.678 

1.2377 

153.689 

1.1128 

123.63 

15 

0.334 

23.162 

0.3079 

18.947 

0.2824 

15.143 

25 

0.1864 

8.  530 

0.1728 

6.969 

0.1594 

5.562 

35 

0.1289 

4.395 

0.1198 

3.588 

0.1109 

2.862 

45 

0.0984 

2.673 

0.0917 

2.1822 

0.0850 

1.740 

55 

0.0796 

1.796 

0.0742 

1.465 

0.0688 

1.168 

80% 

N2 

5 

0.9923 

96.580 

15 

0.2575 

11.754 

25 

0.1464 

4.3113 

35 

0.1021 

2.2175 

45 

0.0783 

1.3475 

55 

0.0636 

0.9046 

' 
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FIG.  21  A 

COMPRESSIBILITY  FACTORS  OF  N2  -  CH4  SYSTEM  AT  32  °F 
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APPENDIX  B 

Determination  of  the  Composition  of  the  Hydrate 
Phase  by  Material  Balance 

As  the  technique  employed  for  the  analysis  of  the 
hydrate  phase  was  likely  to  give  inaccurate  results  due  to 
the  difficulties  in  maintaining  truly  isobaric  conditions  and 
the  possibility  of  the  gas  occlusion  in  the  hydrate  phase,  it 
was  decided  to  develop  a  formula  based  on  the  material  balance 
of  the  system,  for  the  determination  of  the  hydrate  phase 
composition. 

As  stated  earlier,  a  certain  charge  consisting  of 
nitrogen  and  methane  was  compressed  to  a  predetermined  pressure 
by  injecting  mercury  into  the  cell.  During  the  process  of 
hydrate  formation  the  level  of  mercury  in  the  cell  was  main¬ 
tained,  so  as  to  keep  the  volume  of  the  cell  above  the  level 
at  a  constant  value.  As  a  result  of  the  hydrate  formation  the 
pressure  of  the  gas  dropped  from  the  initial  value  of  P1  to 
P  The  samples  of  the  gas  at  the  beginning  and  end  were  analysed 

chromatographically . 
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If  Y±  and  Y2  are  the  mole  fractions  of  nitrogen  at 
the  initial  and  final  conditions;  and  Z^,  Z2  are  the  compressi¬ 
bility  factors  at  the  initial  and  final  conditions,  the  number 
of  moles  of  nitrogen  transferred  from  the  gas  phase  to  the 
hydrate  phase  is, 

M  =  Pl  V°  .  Y  -  .  y2  ( 

Z;lRTo  z2rt0  1 


where  Vq  and  are  the  constant  volume  and  temperature 
of  the  system. 

The  total  number  of  moles  of  gas  transferred  to  the 
hydrate  phase  is. 


Pl  Vo  _  P2  Vo 

ZlRT0  RZ2Tq 


(II) 


Therefore,  the  mole  fraction  of  nitrogen  in  the 
hydrate  phase  is  given  by, 


x  = 


Pl 

Zl 


p2  _ 
1  -  XT 


Z2  *  Y2 


P2 

z2 


(III) 


It  is  assumed  that  the  change  in  volume  of  the  system 
due  to  the  formation  of  hydrate  is  negligibly  small,  and  as 
a  result,  the  volume  of  the  system  remains  constant.  As  the 
solubility  of  the  gas  in  the  water  phase  is  small,  no  correction 
was  made  in  the  derivation  of  the  relation  (III) . 
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The  accurate  estimation  of  'x'  requires  the 
compressibility  factor  data  with  a  high  degree  of  reliability,, 
As  no  adequate  experimental  data  on  the  compressibility  factor 
of  nitrogen-methane  mixtures  were  available,  an  equation  of 
state  (9)  was  employed  to  evaluate  them.  Consequently,  the 
reliability  of  the  calculated  values  of  'x'  using  those  values 
of  compressibility  factors  is  reduced.  The  relation  (III)  was 
therefore  used  with  caution  in  this  investigation. 
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SAMPLE  CALCULATION 

The  application  of  the  Equation  III  is  illustrated 
below  for  the  evaluation  of  the  hydrate  phase  composition.  The 
following  set  of  data  were  obtained  from  a  run  made  at  32°F. 

Initial  pressure  =  p^  =  1223  psia 

Final  pressure  =  P2  =  948  psia 

Mole  percent  ^  in  the  initial  charge  =  Y,  =  56.0 

Mole  percent  ^  in  the  gas  phase  =  Y2  ~  64.3 

at  equilibrium 

Mole  percent  Nj  in  hydrate  phase  =  xact  ~  35.3 

Compressibility  factor  of  the  initial  charge  (Fig.  21) 

=  Z]_  =  0.925 

Compressibility  factor  of  the  gas  at  equilibrium  =  Z2  =  0.955 

-1212.  x  0.439  -  -ZQQ-.  x  0.526 

0.914  0.938 

x  ,  =  =  0.298 

cal . 

1030  _  780 

0.914  0.938 

Mole  percent  N2  in  the  hydrate  phase  by  calculation 
is  29.8 

The  further  application  of  the  Equation  III  is  illust¬ 
rated  through  the  following  table  where  the  calculated  values  of  the 


'  .  '•  - 

“ . -  •  — 
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hydrate  composition  are  compared  with  the  actual  values  of  the 
hydrate  phase  compositions  from  the  experimental  data  at  32°F. 


P1 

psia 

P2 

psia 

Yi 

Mole  %  N2 

Y2 

Mole  %  N2 

xcal. 
Mole  %  N2 

xact. 
Mole  %  N2 

1030.0 

780.0 

49.9 

53,0 

19.0 

20.0 

2240.0 

1860.0 

85.3 

90.0 

62.2 

76.5 

930.0 

524.0 

9.6 

16.0 

2.5 

6.5 

2421.0 

1930.0 

90.6 

92.5 

83.5 

81,0 

1224.0 

948.0 

56.0 

64.5 

29.8 

35.0 

1796.0 

1540.0 

75.9 

81.5 

59.5 

62.0 

2638.0 

2128.0 

91.5 

94.0 

81.8 

86.0 
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APPENDIX  C 


Chromatographic  Analysis 


The  technique  of  gas- liquid  partition  chromatography 
was  employed  for  the  analysis  of  the  gas  samples  consisting  of 
nitrogen  and  methane.  The  Kromo-tog  model  K-2  manufactured  by 
the  Burrel  Corporation  was  used  with  a  2 h  meter  Linde  Molecular 
Sieve  (13-x)  Column.  One  meter  long  column  could  not  separate 
the  oxygen  peak  from  the  nitrogen  peak.  It  was  necessary  to 
have  both  the  peaks  distinctly  separate  so  as  to  know  if  there 
is  any  impurity  in  the  gas  sample.  Helium  gas  with  a  flow 
rate  of  45  ml/min  was  employed  as  a  carrier  gas.  A  cell 
current  of  160  mA  was  used.  Peak  heights  were  used  for  call- 
bration. 


Fluctuations  in  the  room  temperature  had  a  marked 
influence  on  the  peak  heights.  At  times  the  analysis  showed 
a  total  percent  of  more  than  one  hundred.  The  values  of  the 
mole  percents  of  each  of  the  components  were  corrected  to  obtain 

a  total  value  of  one  hundred.  A  similar  correction  was  applied 

.  . .  .  (19,20) 

by  earlier  workers 
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APPENDIX  D 


TABULATIONS 
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TABLE _ I 

EXPERIMENTAL  DATA  ON  INITIAL  HYDRATE -FORMING 
CONDITIONS  FOR  THE  N2-CH4-H2O  SYSTEM 


Composition  -  Mole  % 
Methane 


Pressure 

psia 


Temperature 


100  385,0  32.0 

605,0  40.0 

810,0  45.0 

1257,0  52.7 

1690,0  57.4 

2039,0  60.4 

2916,0  65,4 

4143.0  70 . 1 

87.2  1073.0  49.3 

1350,0  52.5 

2106  c  0  58.2 

2481 „ 0  61.4 

2537.0  63.0 

2831.8  64.1 

2900.0  65.0 

3326.8  67,5 

3576.5  68.5 

4541.4  71.7 

73.1  565.0  32.0 

1298.0  50.3 

1918.0  56.6 

2835.0  62.2 

3769.0  66.5 

4980 . 0  70.3 

50.25  720.0  32.0 

890,0  39.3 

1127.0  43.7 

1521.0  48.5 

2596.0  57.5 

3625.0  62.0 

4814,0  65.6 
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TABLE  I  (Con’t.) 


Composition  -  Mole  % 
Methane 


Pressure 

psia 


Temperature 

°F 


27.2  1155.0  32.0 

1474.0  39.0 
1833.0  44.3 
2472.0  49.5 
2538.0  50.0 
3005.0  53.5 
3647.0  56.5 
4132.0  58.7 

24.0  1250.0  32.0 

1327.0  34.5 
1880.0  41.2 
2530.0  48.0 
3530.0  53.5 
4640.0  58.0 
5215.0  60.7 

10.8  1820.0  32.0 

2300.0  39.2 
2812.0  42.7 
3266.0  46.0 
3745.0  48.0 
4175.0  50.0 

100%  N2  2360.0  32.0 

2485.0  33.0 
2775.0  35.2 
3436.0  38.0 
3655.0  39  o  7 
4150.0  41 o 7 
4390.0  43.0 
5100.0  46.3 
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TABLE  IA 


ISOTHERMAL  DATA  ON  INITIAL  HYDRATE  FORMING  CONDITIONS 


T  = 

Pressure 

psia 

32°F 

y  %  Np 

T  = 

Pressure 

psia 

38°F 

y  %  N2 

T  = 

Pressure 

psia 

44°F 

y  %  N2 

2360.0 

100*0% 

3320,0 

100*0% 

5350.0 

i 

100.0% 

1820.0 

89.2% 

2210.0 

89  *  2% 

3000.0 

89.2% 

1250.0 

76*0% 

1580  .  0 

76.0% 

2850.0 

76.0% 

1155.0 

72.8% 

1400,0 

72.8% 

1830.0 

72.8% 

720.0 

49.8% 

870,0 

VP 

0^ 

CO 

0 

1150.0 

49.8% 

565.0 

26.9% 

680.0 

26.9% 

880.0 

26.9% 

382.0 

0.0% 

520.0 

0.0% 

740.0 

0.0% 

T  =  50°F 

Pressure  y  %  N2 

psia 

T  = 

Pressure 

psia 

56 . 0°F 
y  %  N2 

4200.0 

89.2% 

4100 , 0 

76.0% 

2850.0 

76.0% 

3530 . 0 

72.8% 

2500.0 

72.8% 

2330.0 

49.8% 

1600.0 

49.8% 

1820.0 

26.9% 

1225.0 

26.9% 

1520.0 

0.0% 

1060 . 0 

0.0% 
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TABLE  IB 


ISOBARIC  DATA  ON  INITIAL  HYDRATE  FORMING 
CONDITIONS  FOR  THE  N2~CH4~H20  SYSTEM 


Pressure 

psia 

1000 

1500 

2000 


Temperature 

°F 


40.5 
45.0 

47.5 
49.0 

36.5 

39.5 

48.5 

52.5 
54.0 

35.0 

43.0 

46.0 

53.5 
57.0 
59.0 

33.0 

40.5 

47.5 
50.0 
57.0 

60.5 
63.3 


Mole  %  CH4 


50.3 

73.1 

87.2 

100.0 

24.0 

27.2 

50.3 

73.1 

87.2 

10.8 

24.0 

27.2 

50.3 

87.2 

100.0 

0.0 

10.8 

24.0 

27.2 

50.3 

87.3 

100.0 


2500 
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TABLE  II 


EXPERIMENTAL  PRESSURE  COMPOSITION  DATA 
FOR  THE  N2-CH4~H20  SYSTEM  AT  32°F 


Pressure 

atm. 

Mole  %  N2  in  the 
gas  phase 

Mole  %  N2  in  the 
hydrate  phase 

26.1 

0.0 

0 . 0 

35.7 

16.0 

6.5 

42.5 

31.0 

9.8 

52.8 

53.0 

20 . 0 

64.6 

64.5 

35.0 

76.5 

72.5 

42.5 

105.0 

81.5 

62.0 

115.0 

88.0 

71.0 

126.0 

90.0 

76.5 

131.5 

92.5 

81.0 

144.0 

94.0 

86.0 

160.0 

100.0 

100.0 
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TABLE  III 


EXPERIMENTAL  PRESSURE  COMPOSITION  DATA 
OF  THE  N2-CH4-H20  SYSTEM  AT  44°F 


Pressure 

atm* 

Mole  %  N2  in  the 
gas  phase 

Mole  %  N2  in  the 
hydrate  phase 

50.7 

0.0 

o 

9 

O 

70.5 

35.0 

9.1 

82.6 

46.0 

22.4 

153.5 

75.0 

55.0 

204.0 

CO 

O 

o 

c 

CO 

vO 

249.0 

91.4 

80 . 2 

320.0 

100.0 

100.0 
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TABLE  IV 


EXPERIMENTAL 

PRESSURE  COMPOSITION 

DATA 

FOR  THE 

Nr 

-CH4~H20  SYSTEM 

AT  39 

.  7°F 

Pressure 

atm. 

Mole  %  N2  in 
gas  phase 

the 

Mole  %  N2 
hydrate 

38.1 

0.0 

o 

® 

o 

51.3 

44.0 

00 

0 

o 

00 

o 

G 

O 

63.0 

31  o  0 

102.0 

74.0 

47.0 

119.0 

78.0 

56.  C 

248.5 

100.0 

100.0 

100  - 


TABLE  V 


PRESSURE  COMPOSITION  DATA  FOR  THE  N2-CH4-H20  SYSTEM 

AT  32°F  FROM  THEORY 


y  -  mole  per  cent  nitrogen  in  the  gas  phase 


x  =  mole  per  cent  nitrogen  in  the  hydrate  phase 


P  =  theoretically  predicted  pressure  -  atm. 

Pp  -  corrected  pressure  =  p  (  y  *  UiQ.  +  (1-y)  — — ■ 
Corr  ^  ^  v  *  90  J  18.80 


k  -  equilibrium  constant  =  y/x 


y 

P 

PCorr 

*N2 

*ch4 

0 

18.8 

26.1 

1.000 

5 

19.5 

27.1 

4.460 

01960 

20 

22.3 

32.1 

3.841 

0.845 

40 

27.3 

42.0 

3.140 

0.687 

60 

35.4 

57.0 

2.430 

0.531 

80 

50.1 

85.5 

1.785 

0.363 

100 

90.0 

160.0 

1.000 
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TABLE  VI 


PRESSURE  COMPOSITION  DATA  FOR  THE  CH4-C2H6-H20  SYSTEM 

AT  32°F  FROM  THEORY 


y  «  mole  per  cent  CH4  in  the  gas  phase 

x  =  mole  per  cent  CH4  in  the  hydrate  phase 

P  =  theoretically  predicted  pressure  -  atm. 

,  /  26.1,  / 1  \  5.20  \ 

p_  =  corrected  pressure  =  p  ;  y  *  — — —  +  (1-y;  — ) 

Corr  *  r  \  i  18.8  J  1.10 


k  =  equilibrium  constant  =  y/x 


Y 

X 

P 

p 

Corr 

kCH4 

0 

0.000 

1.10 

5.20 

5 

0.084 

1.13 

5.70 

59.00 

20 

1.450 

1.43 

5.75 

13.80 

50 

5.590 

2.08 

6.30 

8.95 

80 

18.800 

4.26 

9.20 

4.26 

95 

52.750 

11.04 

16.92 

1.80 

100 

100.000 

18.80 

26.10 

1.00 
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TABLE  VII 


PRESSURE  COMPOSITION  DATA  FOR  THE  CH4-C2H4-H20  SYSTEM 

AT  32°F  FROM  THEORY 


y  =  mole  per  cent  methane  in  the  gas  phase 

x  =  mole  per  cent  methane  in  the  hydrate  phase 

P  =  theoretically  predicted  pressure  -  atm. 

pCorr  =  corrected  pressure  =  p(  y  *  (  )  +  (1-y) 

18 . 80 

k  =  equilibrium  constant  =  y/x 


y 

X 

P 

p 

Corr 

kCH4 

0 

0 

0.50 

5.44 

5 

0.179 

0.53 

5.46 

28.80 

20 

0.824 

0.62 

5.58 

24.30 

35 

1.740 

0.76 

5.71 

21.10 

50 

3.170 

0.97 

5.95 

15.80 

65 

5.680 

1.35 

6.35 

11.45 

80 

11.100 

2.22 

7.16 

7.21 

95 

34.400 

6.37 

11.76 

2.76 

100 

100. ooo 

18.80 

26.10 

1.00 

5.44  \ 
0.50 
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TABLE  VIII 


COMPRESSIBILITY  FACTORS  OF  NITROGEN  AND  METHANE 
ALONG  THE  HYDRATE- EQUILIBRIUM  CURVES 


NITROGEN 


Pressure  Temperature  Compressibility  Factor 

atm.  °C 


168.5 

0.5 

1.022 

185.0 

1.5 

1.037 

207.0 

2.5 

1.054 

230.4 

3.5 

1.078 

255.0 

4.5 

1.125 

279.0 

5.5 

1.158 

METHANE 


Pressure  Temperature  Compressibility  Factor 

atm.  °C 


29.6 

1.5 

0.935 

36.1 

3.5 

0.929 

44.0 

5.5 

0.923 

61.3 

8.5 

0.911 

74.9 

10.5 

0.901 

93.5 

12.5 

0.887 

134.5 

15.5 

0.860 

-  io tf  - 


APPENDIX  E 


The  Computer  Programme  for  the  Determination  of 

the  Pressure  Composition  Diagram 


10 


Read,  CMA ,  CMB ,  CNA ,  CNB 
DELMU  =  EXPF  (7.04) 

X  =  0.05 
DELTX  =  0.05 


20  ALPHA  = 

BETA  = 
CHI 

GAMMA  = 

PHI 

ZETA  = 
G 
H 
HA 

C 

D 

If  (C  -  0.0) 


CMA  *  X  +  CNA  *  (1.  -  X) 

CMB  *  X  +  CNB  *  (1.  -  X) 

(3.  *  BETA  +  ALPHA) /( 4.  *  (ALPHA  *  BETA  **  3.)) 

(BETA  **3+3.  *  ALPHA  *  BETA  **  2) 

(4.  *  (ALPHA  *  BETA  **  3)) 

(3.  *  BETA  **  2  +  3  *  ALPHA  *  BETA) 

(6.  *  (ALPHA  *  BETA  **  3)) 

(1.  -DELMU) / (ALPHA  *  BETA  **  3) 

GAMMA  **  2  -  PHI 

PHI  **  3  +  CHI  **  2  -  2  *  PHI *GAMMA  *CHI  +  ZETA*G 
( 4 . *GAMMA *CHI ) / ( 3 .  -  PHI**2  -  ZETA/3 

H  +  SQRTF  (H**2  +  HA**3) 

H  -  SQRTF  (H**2  +  HA**3) 

23,  24,  24 


23  A  =  -  ((-0  **0.3333/2.) 


Go  to  27 


-  io5~  - 


24  A  ( (C) **0.3333/2.) 

27  If  (D  -  O.CQ_  25,  26,  26 

25  B  =  -  ( ( -D) **0.3333/2.) 

Go' to  28 


26  B  ( (D) **0.3333/2.) 
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FIB  =  A  +  B 

U  =  G  +  HB 

V  =  2 .  *G  -  HB 

W  =  4. *U**2  +  3.*ha  -  12 . *G*HB 

PUNCH,  U,V,W 


u 

=  SQRTF 

( U**2) 

V 

=  SQRTF 

( V**2) 

w 

=  SQRTF 

(W**2) 

R1 

=  -  GAMMA 

r2 

=  -  GAMMA 

r3 

=  -  GAMMA 

r4 

=  -  GAMMA 

PUNCH, 

Rl'  R2'  r3 

If  (X  -  1.)  45,  10, 


+  SQRTF  (U)  +  SQRTF 

+  SQRTF  (U)  -  SQRTF 

-  SQRTF  (U)  +  SQRTF 

-  SQRTF  (U)  -  SQRTF 

r4 


( V-H-  SQRTF  (W)  ) 

(V  +  SQRTF  (W)  ) 

(ABSF  (V  -  SQRTF  (W)  )  ) 
(ABSF (V  -  SQRTF (W) ) ) 


45  X  =  X  +  DELTX 

Go  to  20 


CmAj  CriO 
C  N  A,  CNI3 


END 

o.'Tt  the 

owe  ^u*.  L<rt^n  tjm  u-*  A 

Z3  5 


ol  PI  in  a.-nc\  G>^>  Ly 

Cu-r^o  ^cXs?  Ai  /  >1  A.  <j  e  Ca^ji  l  . 

yC  -  lYloU  btn*  r*1  /  rj  ^<5-0  Uj^o  e  _  . 


